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Abstract 

High energy cosmic ray jets from nearby mergers or accretion induced 
collapse (AIC) of neutron stars (NS) that hit the atmosphere can pro- 
duce lethal fluxes of atmospheric muons at ground level, underground 
and underwater, destroy the ozone layer and radioactivate the environ- 
ment. They could have caused most of the massive life extinctions on 
planet Earth in the past 600 My. Biological mutations due to ionizing 
radiations could have caused the fast appearance of new species after 
the massive extinctions. An early warning of future extinctions due to 
NS mergers may be obtained by identifying, mapping and timing all the 
nearby binary neutron stars systems. A warning of an approaching cos- 
mic ray burst from a nearby NS merger/ AIC may be provided by a very 
intense gamma ray burst which preceeds it. 

1 Introduction 

The early history of single-celled organisms during the Precambrian, 4560 to 
570 million years (My) ago, is poorly known. Since the end of the Precambrian 
the diversity of both marine and continental life increased exponentially. Anal- 
ysis of fossil record of microbes, algae, fungi, protists, plants and animals shows 
that this diversification was interrupted by five major mass extinction "events" 
and some smaller extinction peaks [1] . The "big five" mass extinctions oc- 
curred in the Late Ordovician, Late Devonian Late Permian, Late Triassic and 
end-Cretaceous and included both marine and continental life. The largest ex- 
tinction occurred about 251 My ago at the end of the Permian period. The 
global species extinction ranged then between 80% to 95%, much more than, 
for instance, the end-Ordovician extinction 439 My ago which eliminated 57% 
of marine genera, or the Cretaceous- Tertiary extinction 64 My ago which killed 
the dinosaurs and claimed 47% of existing genera [2] . In spite of intensive studies 
it is still not known what caused the mass extinctions, how quick were they and 
whether they were subject to regional variations. Many extinction mechanisms 
have been proposed but no single mechanism seems to provide a satisfactory 
explanation of both the marine and continental extinction levels, the biolog- 
ical extinction patterns and the repetition rate of the mass extinctions [1,2]. 
These include astrophysical extinction mechanisms, such as meteoritic impact 
that explains the iridiun anomaly which was found at the Cretaceous/Tertiary 
boundary [3] but has not been found in all the other extinctions [4], supernova 
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explosions [5] and gamma ray bursts [6], which do not occur close enough at 
sufficiently high rate to explain the observed rate of mass extinctions. 

In this paper we propose that high energy cosmic ray jets (CRJs) from merg- 
ers or accretion induced collapse (AIC) of neutron stars (NS) that hit the atmo- 
sphere of an Earth-like planet can produce lethal fluxes of atmospheric muons 
at ground level, underground and underwater, destroy the ozone layer and ra- 
dioactivate the environment. Nearby NS mergers/ AIC can explain the massive 
extinction on the ground, underground and underwater and the higher sur- 
vival levels of radiation resistant species and terrain sheltered species in the five 
"great" extinctions in the past 600 My. More distant galactic mergers/ AIC can 
cause smaller extinctions. Biological mutations due to ionizing radiations may 
explain the fast appearance of new species after massive extinctions. Intense 
cosmic ray bursts enrich rock layers with detectable traces of cosmogenically 
produced radioactive nucleides such as 129 I, 146 Sm, 205 Pb and 244 Pu. Tracks of 
high energy particles in rock layers on Earth and on the moon may also contain 
records of intense cosmic irradiations. An early warning of future extinctions 
due to neutron star mergers can be obtained by identifying, mapping and timing 
all the nearby binary neutron stars systems. A final warning of an approaching 
CRJ from a nearby neutron stars merger is provided by a very intense gamma 
ray burst a few days before the arrival of the CRJ. 

2 Cosmic Ray Jets From NS Mergers 

Three NS-NS binaries are presently known within the galactic disc; BI9I3+16 [7] 
at a distance of D ~ 7.3 kpc, B203+46 [8] at D ~ 2.3 kpc, and B1534+12 [9] at 
D ~ 0.5 kpc, and B2127+11C [10] in the globular cluster M15 at D ~ 10.6 kpc. 
Although unseen as radio pulsars, the companion stars in these systems have 
been identified as neutron stars from their mass which is inferred via measure- 
ments of the relativistic periastron advance using standard pulse timing tech- 
niques [11]. The continuous energy loss through gravitational radiation brings 
the two NS closer and closer until they merge. It is believed that their final 
merger releases an enormous amount of gravitational binding energy, ~ M Q c 2 , 
in a few ms, in the form of gravitational waves, neutrinos and kinetic energy of 
relativistic ejecta. It was suggested that the neutrinos and/or the relativistic 
ejecta from NS merger/AIC in distant galaxies produce the mysterious gamma 
ray bursts (GRBs), that occur at a rate of about one per day, although the 
mechanism which converts their energy into gamma rays is still not clear [12]. 

It is also believed that due to strong gravitational tidal forces the final NS-NS 
or NS-BH merger proceeds through the formation of an accretion disc. Obser- 
vations seem to indicate that highly collimated jets are ejected by all systems 
where matter is undergoing disc accretion onto a compact central object [13,14]. 
They also indicate that the jet kinetic energy is a considerable fraction of the 
accretion power and that the jets reach enormous distances without being de- 
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fleeted by galactic or intergalactic magnetic fields: The highly relativistic and 
highly collimated jets from active galactic nuclei (AGN), which are believed to 
be powered by mass accretion onto a massive black hole at a typical rate of 
(~ Mq y^ 1 ), reach a distance up to a million light years before disruption [13]. 
Highly collimated relativistic matter that is ejected sporadicly by microquasars 
(superluminal galactic sources such as GRS 1915+105 [14, 15], GRO J165-40 
[16]) and by X-ray sources such as SS 433 [17] and Cygnus X-3 [18], which are 
are close binary systems where mass is accreted onto a neutron star or a stellar 
black hole from a companion star, seem to reach hundreds of light years before 
disruption. The exact mechanism by which the gravitational and electromag- 
netic fields around these accreting and rotating compact objects produce the 
highly relativistic jets is still unknown. However, the instantaneous accretion 
rates (s> M Q s _1 ) and the strength of the magnetic fields which are involved 
in the final stage of merger/AIC of compact stellar objects in compact binary 
systems are probably many orders of magnitude larger than those encountered 
in AGN. Therefore, it is natural to expect that highly relativistic jets are also 
ejected in mergers/ AIC of compact stellar objects, are highly collimated and 
reach hundreds of light years, or more, before disruption. These jets may pro- 
duce the cosmological gamma ray bursts by internal radiation and/or interaction 
with the external medium. After disruption they are isotropized by the galactic 
magnetic field and form the galactic cosmic rays. If they hit an Earth-like planet 
before disruption, they can devastate all forms of life on it. 

3 Constraints From GRBs and Cosmic Rays 

Jets ejected by NS mergers/ AIC in distant galaxies may explain cosmologi- 
cal GRBs [19]. In view of the uncertainties in modeling jet ejection in NS 
merger/AIC, rather than relying on numerical simulations we inferred [19] from 
the observed properties of GRBs that the ejected jets have typical Lorentz fac- 
tors of T ~ 10 3 , beaming angles A£l < 1/100 similar to those observed/estimated 
for AGN and microquasars, and ejected mass AM <~ (dM/<Kl)AVL < 1O _4 M0 
(i.e., released kinetic energy bounded by E K = TAMc 2 < 0.1M Q c 2 ~ 2 x 
10 53 er ff ). 

The finite life times of close binaries due to gravitational radiation emission have 
been estimated [12] from the observed binary period, orbital eccentricity and 
the masses of the pulsar and its companion. They have been used to estimate 
that the NS-NS merger rate in the Milky Way (MW) disc is [20] R M w ~ 
10 -4 — 10~ 5 Beamed ejection from NS mergers/AIC requires merger/AIC 
rate of compact objects in the entire Universe (at GRBs rcdshifts) of the order 
of 10 5 — 10 6 It is larger by 1/AQ <~ 10 2 than that required by spherical 

explosions. It is, however, in good agreement with the updated observational 
and theoretical estimates of the merger rates of compact objects in the Milky 
Way, which yield values in the range 10 5 — 10 6 j/ _1 mergers per Universe [21], 



3 



instead of initial estimates [20] of <~ 10 3 — 10 4 y^ 1 . Thus, the updated estimates 
of the NS-NS merger rate in the Universe are consistent with the observed rate 
[12] of GRBs (approximately one per day) and with the injection rate of cosmic 
rays in the MW that is required in order to maintain a constant cosmic ray flux 
in the MW: 

The escape rate of cosmic rays from the MW requires an average injection 
rate of Qcr ~ 10 41 erg s -1 , in high energy cosmic rays above GeV, in order 
to maintain a constant energy density of cosmic rays in the MW [22]. This 
injection rate can be supplied by NS-NS mergers if the injected jet energy per 
merger is Ek ~ 10 53 er<? and if the NS-NS merger rate in the disc of the MW is 
Rmw ~10- 4 -10" 5 J," 1 : 

Qcr ~ RmwE k - 3 x 10 40 - 3 x 10 41 erg s _1 . (1) 

4 Attenuation of CRJs 

The highly rclativistic jets from quasars and microquasars do not seem to be 
attenuated efficiently in the interstellar or the intergalactic space. This is unlike 
the non relativistic ejecta in supernova (SN) explosions, which is attenuated by 
Coulomb collisions in the interstellar medium (ISM) over a distance of a few 
pc: Moderately energetic charged particles, other than electrons, lose energy in 
neutral interstellar gas primarily by ionization. The mean rate of energy loss 
(or stopping power) is given by the Bcthe-Bloch formula 

- dE/dx « (4^ 2 n e /T 2 e 4 /m e c 2 )[(l/2)ln(2m e c 2 (3 2 T 2 T max / 7 2 ) - 1 ], (2) 

where Ze is the charge of the energetic particle of mass Mj, velocity f3c and total 
energy E — TMiC 2 , n e is the number of electrons per unit volume in the medium 
in atoms with ionization potential I, and T max = 2m e c 2 /3 2 r 2 /(l + 2-/m e /Mi). 
If the interstellar gas around the SN is ionized (by the initial UV flash), then I 
has to be replaced by e 2 /Ro where Rd = (kT/4Tre 2 n e (Z + l)) 1 / 2 is the Debye 
screening length. Thus, for SN ejecta with Z ~ 1 and v ~ 10000 km s _1 in 
an ionized interstellar medium with a typical density of uh ~ 1 cm -3 , and 
kT ~ 1 eV the stopping distance of the ejecta due to Coulomb interactions 
is x ~ E/2(dE/dx) ~ 6 pc. In fact, the stopping of the ejecta by Coulomb 
interactions is consistent with observations of SN remnants, like SN 1006 [23], 
while the assumption that the ionized interstellar medium is glued to the swept 
up magnetic field seems to be contradicted by some recent observations [24]. 
The range increases with energy like /3 4 until nuclear collisions become the 
dominant loss mechanism. The range of nuclei with T ~ 1000 in a typical 
interstellar density of ur ~ 1 cm -3 , is approximately ~ 10 25 cm, i.e., much 
larger than galactic distances. 

Although the galactic magnetic field, H ~ 3 — 5 x 10 -6 Gauss, results in 
a Larmor radius of r L = [3Ymc/qH ~ 10 15 cm for protons with T = 1000, it 
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does not deflect and disperse significantly jets from NS-NS mergers at distances 
smaller than ~ 1 kpc from the explosion. That can be concluded from the fact 
that accretion jets from forming stars, microquasars and AGN reach distances of 
tens, hundreds, and million light years, respectively, without significant deflec- 
tion or attenuation. Probably, because of their high particle and energy densities 
the jets produce internal magnetic fields which shield them from the interstellar 
magnetic field and allow them to follow almost free balistic trajectories in the 
interstellar medium. 



We assume that the ambient interstellar gas is not swept up with the jet. If it 
were, then the jet would not reach even ~ 10 pc. In SN explosions collective 
modes are invoked as the source of the coupling of the ejecta to the interstellar 
medium, required in order to attenuate the SN debris. As mentioned above, 
binary Coulomb interactions are sufficient to produce the observed coupling in 
SN explosions, and coupling through collective modes is not necessarily present. 
Due to internal magnetic fields the jets are highly collimated, not deflected 
and probably reach distances of D ~ 1 kpc, where the internal energy density 
becomes compareble to the external (magnetic and radiation) energy density. 

Unattenuated jets from NS-NS mergers can be devastating to life on nearby 
planets: At a distance of 1 kpc their duration is 



for typical values of T between 1000 and 100, respectively. The time integrated 
energy flux of the jet at D ~ 1 kpc is, typically, ~ 10 12 TeV cmT 2 . Thus, the 
energy deposition in the atmosphere by the jet is equivalent to the total energy 
deposition of galactic cosmic rays in the atmosphere over ~ 10 7 y. However, the 
typical energy of the cosmic rays in the CRJ is ~ 1 TeV per nuclcon, compared 
with ~ 1 GeV per nucleon for ordinary cosmic ray nuclei. Collisions of such 
particles in the atmosphere generate atmospheric cascades where a significant 
fraction of the CRJ energy is converted into "atmospheric muons" through 
leptonic decay modes of the produced mesons. Most of these muons do not decay 
in the atmosphere because of their high energy, unlike most of the atmospheric 
muons which are produced by ordinary cosmic rays. The average number of 
high energy muons produced by nucleons of primary energy E p , which do not 
decay in the atmosphere and reach sea level with energy > at zenith angle 
9 < 7r/2, is given approximately by [25]: 



< JV M >~ (0.0145^[Tey])(^ p /^)°- 757 (l - E^/E p ) 5 - 25 /cos9 . (4) 



Thus a jet with energy of about 1 TeV per nuclcon at a distance of 1 kpc 
produces at sea level a flux of atmospheric muons of 



5 Mass Extinctions By CRJs 



St ~ D/2cT 2 ~ 1 



day — 2 months 



(3) 



/„(> 3 GeV) ~ 10 12 cm 



-2 



(5) 
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Such muons deposit energy in matter via ionization. Their energy deposition 
rate is [26] -dE/dx > 2 MeV g~ r cm~ x . The whole-body lethal dose from 
penetrating ionizing radiation resulting in 50% mortality of human beings in 30 
days [26] is < 300 rad ~ 2 x 10 10 /(dE/dx) ~ 10 10 cm~ 2 where dE/dx is in rate 
is in MeV <7 _1 cto _1 units. The lethal dosages for other vertebrates can be a few 
times larger while for insects they can be as much as a factor 20 larger. Hence, a 
CRJ at D ~ 1 kpc which is not significantly dispersed by the galactic magnetic 
field produces a highly lethal burst of atmospheric muons. Because of muon 
penetration, the large muon flux is lethal for most species even deep (hundreds 
of meters) underwater and underground, if the cosmic rays arrive from well 
above the horizon. Thus, unlike the other suggested extraterrestrial extinction 
mechanisms, a CRJ which produces a lethal burst of atmospheric muons can 
explain also the massive extinction deep underwater and why extinction is higher 
in shallow waters. 

Although half of the planet is in the shade of the CRJ, planet rotation exposes 
a larger fraction of the planet surface to the CRJ. Additional effects increase 
the lethality of the CRJ over the whole planet. They include: 

(a) The pollution of the environment by radioactive nuclei, produced by 
spallation of atmospheric and surface nuclei by shower particles. Using the 
analytical methods of [27], we estimate that for an Earth-like atmosphere, the 
flux of energetic nucleons which reaches the surface is also considerable, 

I p (> 100 MeV) ~ I n (> 100 MeV) ~ 10 10 cm,- 2 . (6) 

Global winds spread radioactive gases in a relatively short time over the whole 
planet. 

(b) Depletion of stratospheric ozone by the reaction of ozone with nitric ox- 
ide, generated by the cosmic ray produced electrons in the atmosphere (massive 
destruction of stratospheric ozone has been observed during large solar flares 
which produced energetic protons [28]). 

(c) Extensive damage to the food chain by radioactive pollution and massive 
extinction of vegetation and living organisms by ionizing radiations (the lethal 
radiation dosages for trees and plants are slightly higher than those for animals 
but still less than the flux given by eq. 5 for all except the most resilient species). 

6 Signatures of CRJ Extinction 

The biological extinction pattern: The biological extinction pattern due to 
a CRJ depends on the exposure and the vulnerability of the different species 
to the primary and secondary effects of the CRJ. The exposure of the living 
organisms to the muon burst depends on the intensity and duration of the CRJ, 
on its direction relative to the rotation axis of Earth (Earth shadowing), on 
the local sheltering provided by terrain (canyons, mountain shades) and by un- 
derwater and underground habitats, and on the risk sensing/assessment and 
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mobility of the various species. The lethality of the CRJ depends as well on 
the vulnerability of the various living species and vegetation to the primary 
ionizing radiation, to the drastic changes in the environment (e.g., radioactive 
pollution and destruction of the ozone layer) and to the massive damage and 
radioactive poisoning of the food chain. Although the exact biological signa- 
ture may be quite complicated, and somewhat obscured in fossil records (due to 
poor or limited sampling, deterioration of the rocks with time and dating and 
interpretation uncertainties because of bioturbational smearing) it may show 
the general pattern expected from a CRJ extinction. Indeed, a first examina- 
tion of the fossil records suggest that there is a clear correlation between the 
extinction pattern of different species, their vulnerability to ionizing radiation 
and the sheltering provided by their habitats and the environment they live in. 
For instance, insects which are less vulnerable to radiation, were extinct only 
in the greatest extinction - the end-Permian extinction 251 My ago. Even then 
only 8 out of 27 orders were extinct compared with a global species extinction 
that ranged between 80% to 95% [4]. Also plants which are less vulnerable 
to ionizing radiation suffered lower level of extinction. Terrain, underground 
and underwater sheltering against a complete extinction on land and in deep 
waters may explain why certain families on land and in deep waters were not 
extinct even in the great extinctions, while most of the species in shallow wa- 
ters and on the surface were extinct [4]. Mountain shadowing, canyons, caves, 
underground habitats, deep underwater habitats and high mobility may also 
explain why many species like crocodiles, turtles, frogs, (and most freshwater 
vertebrates), snakes, deep sea organisms and birds were little affected in the 
Cretaceous/Tertiary (K/T) boundary extinction which claimed the life of the 
big dinosaurs and pterosaurs. In particular, fresh underground waters in rivers 
and lakes are less polluted with radioisotopes and poisons produced by the CRJ 
than sea waters and may explain the survival of freshwater amphibians. 
Geological Signatures: The terrestrial deposition of the primary CRJ nu- 
clides or the production of stable nuclides in the atmosphere or in the surface 
by the CRJ is too small to be detectable. In particular, the proposed mech- 
anism cannot explain the surface enrichment at the K/T boundary by about 
3 x 10 5 tons of iridium. Alvarez et al. [3] suggested that the impact of ex- 
traterrestrial asteroid, with an Ir abundance similar to that observed in early 
solar system Chondritic meteorites whose Ir abundance is larger than crustal Ir 
abundance by ~ 10 4 , could cause the Ir anomaly and explain the mass extinc- 
tion at the K/T boundary. However, no significant Ir enrichment was found in 
all other mass extinctions. Moreover, other isotopic anomalies due to meteoritic 
origin have not been found around the K/T boundary; in particular the As/Ir 
and Sb/Ir ratios are three orders of magnitude greater than chondritic values 
but are in accord with a mantle origin [29]. Extensive iridium measurements 
showed that the anomaly does not appear as a single spike in the record, indica- 
tive of an instantaneous event, but rather occur over a measurable time interval 
of 10 to 100 ky or possibly longer [30]. That and the high abundance of Ir 
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in eruptive magma led to the suggestion that the iridium anomaly is due to a 
global volcanic activity over 10 to 100 kyr at the K/T boundary [30,31] which 
also caused the K/T extinction: Eruptions have a variety of short term effects, 
including cooling from both dust and sulfates ejected into the stratosphere, acid 
rain, wildfires, release of poisonous elements and increase in ultraviolet radia- 
tion from ozone-layer depletion. But, examination of major volcanic eruptions 
in the past 100 My have shown that none of them greatly affected the diversity 
of regional and global life on land or in the oceans [4] . 

A CRJ will enhance the abundance of stable cosmogenic isotopes in the 
geological layer corresponding to the CRJ event, but, the enrichment may be 
negligible compared to their accumulation through long terrestrial exposure of 
the geological layers to galactic cosmic rays prior to the CRJ. However, CRJ 
enrichment of sediments with unstable radioisotopes of mean lifetimes much 
shorter than the age of the solar system, t <C t & w 4570My, but comparable to 
the extinction times, may be detectable through low traces mass spectrometry. 
In particular, fission of long lived terrestrial nuclei, such as 238 U and 232 Th, by 
shower particles, and capture of shower particles by such nuclei, may lead to 
terrestrial production of, e.g., 129 I with r = 15 My, 146 Sm with r = 146 My, 
205 Pb with t = 43 My and 244 Pu with r = 118 My, respectively. These 
radioisotopes may have been buried in underwater sediments and underground 
rocks which were protected from further exposure to cosmic rays. The main 
background to such a CRJ signature is the continuous deposition by cosmic 
rays and by meteoritic impacts on land and sea. Cosmic rays may include 
these trace radioisotopes due to nearby sources (e.g., supernova explosions) and 
because of spallation of stable cosmic ray nuclei in collisions with interstellar 
gas. Meteorites may include these trace elements due to a long exposure in 
space to cosmic rays. 

Finally, large enhancement of TeV cosmic ray tracks in magma from volcanic 
eruptions coincident with extinctions may also provide fingerprints for CRJ 
extinctions. 

7 Rate of Mass Extinctions 

Assuming that the spatial distribution of NS binaries and NS mergers in the 
MW follow the distribution of single pulsars [] , 

dN <xe- R2 / 2R °e-W h RdRdz, (7) 

with a disc scale length, Ro ~ 4.8 kpc, and a scale height, h > 0.5 kpc perpendic- 
ular to the disc and independent of disc position, we find that the average rate 
of CRJs from NS-NS mergers that reach planet Earth from distances < 1 kpc, is 
~ 10~ 8 y^ 1 . It is consistent with the 5 big extinctions which have occurred dur- 
ing the last 600 My in the Paleozoic and Mesozoic eras. The relative strengths 
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of these extinctions may reflect mainly different distances from the CRJs. Be- 
yond ~ 1 kpc from the explosion the galactic magnetic field begins to disperse 
the CRJ and suppresses its lethality. Such CRJs, if not too far, can still cause 
partial extinctions at a higher rate and induce biological mutations which may 
lead to the appearance of new species. 

The galactic rate of SN explosions is <~ 100 y^ 1 . The range of debris from 
SN explosions in the interstellar medium is shorter than 10 pc. The rate of 
SN explosions within a distance of 10 pc from Earth which follows from eq. 7 
is Rsn(< 10 pc) m 10~ 10 y~ 1 . High energy cosmic rays which, perhaps, are 
produced in the SN remnant by shock acceleration, carry only a small fraction 
of the total explosion energy and arrive spread in time due to their diffusive 
propagation in the interstellar magnetic field. Also neutrino and light emissions 
similar to those observed in SN1987A, at a distance of a few pc cannot cause a 
mass extinction. 

8 Conclusions 

Cosmic Ray Bursts from neutron star mergers may have caused the massive 
continental and marine life extinctions which interrupted the diversification of 
life on our planet. Their rate is consistent with the observed rate of mass 
extinctions in the past 570 My. They may be able to explain the complicated 
biological and geographical extinction patterns. Biological mutations induced 
by the ionizing radiations which are produced by the CRJs may explain the 
appearance of completely new species after extinctions. A first examination 
suggests a significant correlation between the biological extinction pattern of 
different species and their exposure and vulnerability to the ionizing radiation 
produced by a CRJ. The iridium enrichment around the Cretaceous/Tertiary 
extinction that claimed the life of the dinosaurs and pterosaurs cannot be due 
to a CRJ. It may have been caused by intense volcanic eruptions around that 
extinction. Isotopic anomaly signatures of CRJ extinctions may be present in 
the geological layers which recorded the extinctions. Elaborate investigations 
of the effects of CRJs from relatively nearby neutron star mergers and their 
biological, radiological and geological fingerprints are needed before reaching a 
firm conclusion whether the massive extinctions during the long history of planet 
Earth were caused by CRJs from neutron stars mergers. If nearby neutron star 
mergers are responsible for mass extinctions, then an early warning of future 
extinctions due to neutron star mergers can be obtained by identifying, mapping 
and timing all the nearby binary neutron stars systems. A final warning for an 
approaching CRJ from a nearby neutron-stars merger will be provided few days 
before its arrival by a gamma ray burst produced by the approaching CRJ. 
Acknowledgement: This research was supported in part by the Technion fund 
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